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Abstract 

Analytical  expression  of  electric  power  was  deduced  in  case  of  the  large-scale  cylindrical  thermoelectric  tubes  exposed  to  two  thermal 
fluids.  The  output  powers  of  the  proposed  six  geometrical  arrangements  were  mathematically  solved  from  heat  transfer  theory,  and  compared 
with  the  flat  panel  systems.  The  maximum  output  was  the  largest  in  the  ideally  cooled  systems.  In  the  other  realistic  systems,  it  was  the  largest 
for  the  system  of  the  counter  flow  with  a  single  thermoelectric  tube.  The  multiplication  of  thermoelectric  panels  can  shorten  significantly 
the  device  length,  although  the  output  from  the  multi-tubes  decreases  only  a  little.  For  example,  the  double  helical  tubes  (T2CH  system) 
can  generate  96%  output  by  35%  length,  compared  with  the  single  tube  in  the  counter  flow  (TIC  system). 

©  2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  a  large-scale  power  plant,  a  thermal  energy  is  carried 
to  the  power  generator.  When  it  is  extracted  by  a  fluid  and 
passed  to  the  thermoelectric  generator,  a  hot  fluid  offers  the 
heat  to  the  hot  junctions  in  the  thermoelectric  modules.  See- 
beck  effect  can  generate  the  thermoelectric  power  depending 
on  a  temperature  difference  between  the  hot  and  cold  junc¬ 
tions.  The  thermoelectric  motive  force  in  the  open  circuit  is 
the  sum  of  multiplication  of  the  relative  Seebeck  coefficient 
and  the  temperature  difference  AT  over  all  the  serial  connec¬ 
tions.  The  problem  to  obtain  the  larger  power  is,  therefore, 
how  to  give  the  larger  AT  to  all  the  thermoelectric  panels 
existing  between  two  hot  and  cold  fluids  [1-4]. 

The  heat  applied  from  a  hot  fluid  is  transferred  to  the 
panel  surface,  into  the  thermoelectric  panel,  and  finally  to  the 
cold  fluid  at  another  surface.  The  two  fluids  are  resultantly 
warmed  or  cooled  along  the  flow  paths,  and  their  tempera¬ 
ture  profiles  through  the  path,  T(x),  change  as  a  function  of 
position,  x.  This  problem  is  well  known  in  heat  exchanger 
through  an  isolator  [5-7].  Additionally  as  the  special  fea¬ 
ture  in  the  thermoelectric  power  generation,  the  heat  transfer 
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due  to  Peltier  effect  and  Joule  heating  should  be  considered 
[1,8,9]. 

The  purpose  of  this  work  is  to  show  the  mathematical 
expression  of  the  electric  power  extractable  from  the  ther¬ 
moelectric  cylindrical  tubes  that  are  heated  by  a  hot  fluid 
and  cooled  by  another  cold  fluid.  Fig.  1  shows  an  example 
of  this  concept.  The  system  of  Fig.  1  consists  of  three  ther¬ 
moelectric  tubes,  in  which  two  hot  fluids  and  the  two  cold 
fluids  flow  as  counter  flow. 

The  previous  studies  showed  that  the  maximum  of  output 
power  exists  corresponding  to  a  certain  module  size  [8-11]. 
It  is  because  the  longer  serial  connections  of  thermoele¬ 
ments  can  generate  the  higher  voltage,  but  the  extremely 
longer  panels  have  the  higher  electrical  resistivity  and  the 
temperature  difference  AT  becomes  smaller.  For  example, 
Kyono  et  al.  evaluated  the  case  when  a  single  thermoelectric 
tube  was  isothermally  cooled  by  the  huge  amount  of  coolant 
[10].  They  reported  that  the  maximum  of  150kW  existed, 
although  this  optimum  length  was  evaluated  as  a  few  hun¬ 
dreds  km  for  FeAlSi-FeAl  couple.  However,  this  optimum 
size  of  the  thermoelectric  generator  can  be  compacter  in  a 
few  tens  metre  by  using  the  high  performance  semiconduc¬ 
tors  such  as  Bi2Te3  [10,11]. 

It  is  expected  that  the  multi-tubes,  as  illustrated  in  Fig.  1, 
can  shorten  the  tube  length  and  make  the  system  more 
compact  and  suitable  for  the  large-scale  power  genera¬ 
tion,  because  the  previous  study  showed  that  the  necessary 
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Thermoelectric  tubes 


Adiabatic  wall 

Fig.  1.  Cylindrical  multi-tube  system. 


panel  area  could  be  minimized  when  the  thermoelectric 
flat  multi-panels  were  connected  three-dimensionally  [11]. 
When  the  heat  loss  through  a  thermoelectric  tube  is  recov¬ 
ered  by  another  fluid,  and  when  the  captured  heat  is  passed 
to  another  thermoelectric  panel,  we  expect  that  the  thermal 
energy  can  be  efficiently  recovered  in  a  restricted  space 
[2-4,8-11].  The  cylindrical  tubes  consisting  of  the  thermo¬ 
electric  modules  are  suitable  for  circulation  of  the  industrial 
fluids,  such  as  exhaust  hot  gas  or  coolant  water. 

The  analytical  method  for  flat  multi-panels  [11]  and  that 
for  a  cylindrical  tube  [10]  are  here  combined,  and  the  heat 
transfer  through  the  thermoelectric  multi-tunes  and  the 
temperature  change  T(x)  along  the  fluid  path  with  cylin¬ 
drical  symmetry  are  analyzed  mathematically  to  discuss 
the  maximum  output  power.  Here  we  propose  six  types  of 
thermoelectric  power  generation  system  with  cylindrical 
multi-tubes. 


fluids  are  isolated  by  the  panel,  and  flow  along  the  both 
cylindrical  surfaces  of  the  panel.  Practically,  as  shown  in 
Fig.  2(a),  there  exist  the  electrodes  connecting  the  thermoele¬ 
ments,  the  insulator  between  the  elements,  the  protective  in¬ 
sulating  film  on  the  electrodes,  and  the  fins  that  enhance  the 
heat  exchange.  Here  for  simplicity  we  neglected  these  ef¬ 
fects  on  heat  transfer,  and  considered  them  only  as  the  heat 
transfer  coefficient  of  the  modules. 

2.2.  Directions  of  fluids 

The  number  of  the  tubes,  n ,  the  directions  of  hot  and  cold 
fluids,  and  the  path  shape  for  fluids  are  classified  systemat¬ 
ically  as  shown  in  Fig.  3. 

Firstly,  three  types  of  fluid  direction  are  considered,  de¬ 
pending  on  the  direction  of  cold  fluid  at  the  outlet  of  the 
corresponding  hot  fluid:  parallel  flow  type  (P),  counter  flow 
type  (C),  and  isothermal  type  (I).  We  assume  only  one  path 
for  each  fluid.  Namely,  we  set  the  boundary  conditions  that 
each  system  has  only  two  paths  connecting  in  series.  No 
branched  paths  are  considered,  although  it  was  one  of  the 
efficient  methods  as  reported  previously  [11]. 

The  parallel  flow  type  is  the  system  that  the  two  fluids 
flow  in  the  same  direction,  between  which  a  thermoelectric 
panel  is  inserted.  The  fluids  flow  in  the  opposite  direction  at 
the  counter  flow  type.  These  two  types  are  often  selected  for 
heat  exchangers  using  fluids  [5-7].  The  isothermal  type  is 
the  extreme  case  that  whole  a  surface  of  the  most  outside  of 
the  system  is  kept  at  a  constant  temperature  by  the  infinitely 
large  heat  bath.  This  case  for  a  single  tube  (Til)  was  partially 
analyzed  by  Kyono  et  al.  [10]. 


2.  Basic  assumptions  and  models 

2.7.  Cylindrical  thermoelectric  tube 

It  is  assumed  that  all  of  our  thermoelectric  modules  are 
cylindrical,  and  that  they  consist  of  the  thermoelements  with 
a  single  layer  of  n-type  p-n  junctions,  as  illustrated  in 
Fig.  2(a).  The  thermoelements  are  homogeneously  aligned 
perpendicular  to  heat  flow,  combined  tightly  without  open 
space,  and  connected  electrically  in  series.  The  hot  and  cold 


Fig.  2.  Simplification  of  thermoelectric  tube,  (a)  Practical  model  and  (b) 
simplified  model. 
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Fig.  3.  Classification  of  thermoelectric  systems. 
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Table  1 


Conditions  for  the  analyzed  systems 


System 

Hydrodynamic  conditions 

Boundary  conditions 

In  Fig.  3 

TIP 

Mi  >  0,  M2  >  0 

Ti  (0)  =  Tf,  72(0)  =  Tf 

Fig.  3(a) 

TIC 

Mi  >  0,  M2  <  0 

Ti  (0)  =  Tf,  r2(Z)  =  Tf 

Fig.  3(b) 

Til 

Mi  =  M  >  0 

Tx  (0)  =  Tf,  T2(x)  =  Tf  for  all  jc 

Fig.  3(c) 

T2CM 

Mi  >  0,  M2  —  —Mi,  M3  =  Mi 

T\  (0)  =  Tf,  Tx  (/)  =  72(/),  ?3  (0)  =  Tf 

Fig.  3(d) 

T2CH 

Mi  >  0,  M2  <  0,  M3  Mi 

TX  (0)  =  Tf,  r2  (0  =  Tf,  73(0)  =  Ti(l) 

Fig.  3(e) 

T2I 

Mi  >  0,  M2  <  0 

Ti  (0)  =  Tf,  r2  (0  =  Tf,  73  (x)  =  Tf  for  all  x 

Fig.  3(f) 

2.3.  Path  shape 

Secondly  the  path  shapes  of  fluids  are  classified  into  me¬ 
andering  type  (M)  and  helical  type  (H),  as  shown  in  Fig.  3. 
The  meandering  type  and  the  helical  type  aim  at  the  heat 
recycling  and  at  the  homogeneous  AT  for  all  the  panels,  re¬ 
spectively.  The  path  connection  is  illustrated  as  broken  lines 
in  Fig.  3(e).  The  seamless  three-dimensional  connection  is 
realized  by  bending  the  tube. 

The  analyzed  tube  (T)  systems  were  named  after  the  num¬ 
ber  of  panel  sheets  ( n ),  the  direction  of  fluid  flow  (P,  C  or  I) 
and  the  path  shape  (M  or  H).  The  hydrodynamic  conditions 
and  the  boundary  conditions  are  listed  in  Table  1 . 


3.  Equations  for  output  power 

3.1.  Deduction  of  equations 

Our  previous  approach  for  multi-panels  is  applied  for  the 
multi-tubes  [11].  The  modification  of  heat  transfer  through 
the  cylindrical  thermoelectric  tube  is  considered  [10].  The 
circumferential  homogeneous  temperature  is  assumed  that 
the  turbulent  flow  is  well  developed  inside  the  long  path. 

When  the  system  consists  of  n  tubes,  a  set  of  n  +  1  simul¬ 
taneous  derivative  equations  can  be  written  from  the  heat 
transfers  through  tubes, 

d7\(x) 

MlCp4^r  =  TKl(Tl(x)  ~  Tl{x))  (1) 


shown  in  Fig.  4.  A  is  the  average  heat  conductivity  of  the 
tube,  given  by, 

^  _  ^p0p  +  ^n0n 
0p  +  0n 

where  Ap,  An,  0p  and  0n  are  the  heat  conductivity  and  the 
angle  of  p-  and  n-type  elements,  respectively.  For  simplicity, 
heat  transfer  by  Peltier  effect  and  exothermal  heat  by  Joule 
effect  are  neglected.  This  assumption  leads  to  an  overesti¬ 
mation  for  output  power,  although  this  overestimation  is  not 
significant,  as  we  will  report  separately. 

The  equation  is  solved  for  x  by  setting  the  boundary  con¬ 
ditions  listed  in  Table  1.  The  temperature  6ij(x)  at  the  sur¬ 
face  of  the  panel  i  facing  to  the  fluid  j  is  then  written  using 
the  solutions  of  simultaneous  derivative  equations,  Tfx). 

0u(x)  =  Ti(x)  -  Ti(x )  -  7}+i(x))  (6) 

hi 

0u+ i(*)  =  Ti+ i(x)  +  Tfx )  -  Ti+i(x))  (7) 

hi+ 1 


The  electromotive  force,  E ,  is  the  summation  of  the  temper¬ 
ature  difference  over  all  the  panels, 


0ij+i(x))dx 


where  a  is  the  difference  between  the  Seebeck  coefficients  of 
p  and  n  elements,  n 0  and  nx  are  the  number  of  thermoelectric 
pairs  in  a  circumferential  circulation,  and  the  number  density 


dTi(x) 

=F  Ki(Ti(x)  -  Ti+i(x))  (1  <  i  <  n)  (2) 
dTn+itx) 

M„+1Cp,„+i  =  ±Kn(Tn(x)  -  Tn+ i(x))  (3) 


where  x  is  the  position  along  the  tube  length.  M;,  Cpj  and  Ti 
are  mass  flow  rate,  heat  capacity  and  temperature  for  fluid  i, 
respectively.  The  choice  of  ±  depends  on  the  path  condition. 
K[  is  the  over-all  heat  transfer  coefficient  through  the  tube  i 
in  the  direction  perpendicular  to  x-axis, 


2tt 


(1/ hijTij)  +  (ln(r;+i/r;)/A)  +  (1/ hi+ijri+ij) 


r  and  h  are  the  radius  of  the  tube  and  the  heat  transfer 
coefficient  between  the  fluid  and  the  tube,  respectively,  as 


Fig.  4.  Dimensions  of  thermoelectric  tubes. 
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of  the  pairs  in  the  direction  x,  respectively.  Note  that  = 
27t/(0p  +  0 n).  The  electric  resistance,  Ri,  is  given  as, 


)  (9) 

\  <Pp  4>nJ  \  rU  ) 

where  pv  and  pn  are  the  specific  resistivity  of  p  and  n  ele¬ 
ments.  The  output  power,  P,  is  optimized  by  balancing  the 
internal  and  external  resistance, 

P  =  4 TURi  (10) 

This  work  uses  hereafter  this  P,  where  the  internal  resistance 
is  equal  to  the  external  resistance. 

3.2.  Output  power  of  six  systems 

The  simultaneous  equations  are  solved  at  the  TIP  system, 
where  one  sheet  of  panel  is  exposed  in  parallel  flow.  The 
output  was  calculated  by  Eq.  (10)  as, 

n0^{M1M2Cp,iCp,2a(rijn  -  ^n)}2 


Ptip  = 


X 


16(rn  +  rv)l 

{e(-(1/MlCP,1  )-^/m2cv,2))ki  _  { j 
^lCp,l  +  ^2<7p,2 


n  2 


(ID 


where  rn  and  rp  are  the  resistivity  for  n-  and  p-type,  respec¬ 
tively,  in  one  unit  of  thermocouple.  The  output  Ptip  had 
been  maximized  using  the  conditions  that, 

27Ty//Ap  pn 

<t>n  = - ,  - .  and 

W0(yAnPp  +  yAp/)n) 

J,  _  2 

0p  —  i -  i -  (12) 

^nPp  T  y hpPn) 

and  the  new  parameter  K%\  was  defined  relating  with  heat 
flow  through  the  panel  as, 

h\\K\r\\  +  h\2K\r\2  -  2jih\\h\2r\\r\2 


K  al  = 


jxh\\h\2K\r\\r\2 


ln(ri2/rn) 

ji\\ 


(13) 


The  same  optimized  conditions  for  0n  and  0p  were  also  used 
in  the  following  calculations. 

The  solutions  for  TIC  system  were  generally  given  as, 

n0^{M1M2Cp,iCp,2«(7]jn  -  Tcin)}2 


Ptic  = 


x 


-i  2 


16(rn  +  rp)/ 

MiCp4e(“(1/MlCP’l)“(1/M2CP’2))^/  +  M2Cp,2 

(14) 


When  MiCpj  =  —  M2Cp? 2  =  MCp,  however,  the  special 
care  was  needed  in  solving  the  differential  equations.  The 
output  power  of  this  case, 

n^nAMCpdiT™  -  T™)}21  r  -  n2 


^TIC  = 


16(rn  +  rp) 


gal 
_MCp  +  Kl_ 


(15) 


was  equal  to  the  infinite  expression  of  Eq.  (13)  when  M \  Cp_  i 
approaches  to  — M2Cp  2. 

The  solutions  for  Til  system  was  given  as, 


nn\i2 


^TII  = 


n4,nx{MCpa(T™  -  Tcm)} 
16(rn  +  rp)l 


[ATal(l  -  e“(A7/'MCP,)l2 

(16) 


Eq.  (16)  is  the  infinite  expression  when  M2Cp52  is  infinitely 
large  in  either  Eq.  (11)  or  (14).  It  is  natural  because  the  Til 
system  holds  the  infinitely  large  heat  bath,  which  is  made 
possible  by  supply  of  an  infinitely  large  amount  of  fluid. 
Eq.  (16)  essentially  agreed  with  the  previous  report  [10], 
although  the  definitions  of  the  symbols  were  different. 

In  the  solutions  for  the  double-tubes  systems,  we  will 
report  only  the  case  of  \M\CV:\  \  =  |M2Cp52l  =  IM3 Cp?3|  = 
MCp  for  simplicity, 


^T2CM  = 


 nfiAMCpaijt  ~  CT 


2(r nl  T  Pn2  T  Ppl  +  P p2)^ 

[sin  h  (X)  {/Cai  yZT  sinh(2Q 
+  ^a2V^2COSh(X)}]2 


X 


VZT cosh(22Q  +  ^[K~2  sinh(2X) 


(17) 


where  K&  and  X  were  defined  as, 

h  22^2r22  +  ^23^2^23  —  277^22^23  r22r23 


^a2  = 


71/722^23  ^2r22r23 


and 


X  = 


In(r23/P22) 

2 


l^KfK~2 

2MCn 


(18) 


(19) 


respectively.  Using  the  same  parameters, 


Et2CH  = 


n,pnx{MCpa(T™  -  Tcin)sinh(X)}2 
4(rni  T  rn2  +  Tpj  T  r*p2)(/fai  Ka2)l 

{(KxKal  +  K2Ka2)  cosh(X) 

+  jKOG(Kal  -  Kal)  sinh(X)} 

x 


_ 


TRSKi  cosh(2Z) 

+  (ATi  +  K2)  sinh(2X) 


(20) 


and 


Et2I  = 


n^nAMCpaiT™  -  T™)V 
16(rni  +  rn2  +  rpi  +  rp2)l 

jKAAKi  +  A'2)A'a2(cosh(Z)  -  e~Y} 

+  (2KiKal  +  K2Ka2)  sinh(Z) 

x 


_ 


V^2(4gi  +  g2)  cosh(Z) 
+  (2 Ki  +  K2)  sinh(Z) 


(21) 


where 


Y  = 


Kl 

Ki 


X  and 


Z  =  J4+—X 
K 1 


(22) 
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The  obtained  expressions  contain  the  common  term  of 
n(f)nx{M\M2Cv,\Cv,2w(T^i  —  Tcm)}2.  The  other  parts  de¬ 
pend  on  the  geometry  of  the  system,  electric  resistivity 
and  heat  transfer  coefficient  K(.  Ki  depends  on  the  radius 
of  thermoelectric  panel  in  the  cylindrical  system.  This  is 
unable  to  simplify  the  expressions,  while  the  output  power 
of  the  flat  multi-panels  systems  could  be  expressed  by  the 
non-dimensional  parameters  [11].  Therefore,  the  obtained 
output  for  six  systems  will  be  numerically  compared  using 
the  thermophysical  values  of  Bi2Te3  semiconductors. 

4.  Physical  properties  and  conditions 

The  thermophysical  properties  of  the  BiTe  thermoelec¬ 
tric  elements  in  the  literature  scattered  due  to  the  impurities, 
their  temperature  dependencies  and  the  difference  in  prepa¬ 
ration  [1,12].  Table  2  shows  a  typical  set  of  properties  for  an 
identical  sample  at  room  temperature  [1],  and  used  for  this 
work.  Table  3  shows  the  fluid  properties  [6,7,13]  and  the 
parameters  for  thermoelectric  tubes.  The  thickness  of  ther¬ 
moelectric  elements  was  set  as  0.05  m  based  on  the  previ¬ 
ous  discussion  [11].  The  most  inner  radius  was  set  as  0.1  m, 
and  the  path  width  (the  gap  between  cylindrical  two  tubes) 
0.1  m.  The  other  thermophysical  properties  and  conditional 

Table  2 


Specific  values  of  thermoelectric  elements  used  here 


Materials 

Seebeck 
coefficient, 
a  (|iV/K) 

Resistivity, 
p  (fx£2m) 

Thermal 
conductivity, 
X  (W/Km) 

Figure  of 
merit,  Z  (1/K) 

Bi-54.3  at.% 

162 

5.55 

2.06 

Te  (p) 

Bi-64.5  at.% 

-240 

10.1 

2.02 

2.605  x  10“3 

Te  (n) 

Table  3 

Parameters  for  thermoelectric  power  generation  system 

Variables 

Values  used  in 

this  work 

Thermoelectric 

Length,  / 

Variable 

device 

Thickness  of  tube,  d 

0.05  m 

Radius  of  inner  tube,  r\ 

0.1  m 

Distance  between  two  tubes 
(thickness  of  fluid  path) 

0.1  m 

Number  density  of  pairs,  nx 

50/m 

Number  of  pairs,  n(j} 

100/cy  finder 

Thermal  sources 

Hot  source 

N2  gas  (inlet 

7]jn  =  500  K) 

Cold  source 

N2  gas  (inlet 

Tcin  =  300  K) 

Gas  properties 

Specific  heat,  Cp^ot  —  Gp  COOi 

1044  J/(kg  K) 

(at  400  K, 

0.1  Pa) 

Mass  flow  rate,  Mhot  =  Mcooi 

1.00  kg/s 

Heat  transfer  coefficient, 

/*hot  —  ^cool 

500  W/(m2  K) 

parameters  are  identical  with  the  previous  work  [11],  and 
the  temperature  dependencies  of  these  values  are  not  con¬ 
sidered. 


5.  Comparison  of  multi-tube  systems 

Fig.  5  shows  the  output  power  thus  evaluated  for  six  sys¬ 
tems.  The  output  power  of  each  system  became  the  maxi¬ 
mum  at  a  certain  length.  This  maximum  output  power,  Pmax, 
and  its  system  length,  Lmax,  are  listed  in  Table  4.  Pmax  was 
the  largest  for  Til  system,  i.e.  P™ =  9.49  kW  at  Lmax  = 
48.5  m.  However,  the  isothermal  heat  sink  in  the  Til  system 
is  an  ideal  case,  and  it  is  practically  impossible  to  supply  the 
fluid  at  the  infinitely  fast  rate  in  order  to  keep  isothermal. 
When  we  feed  the  hot  and  cool  fluids  at  the  same  finite  rate 
(i.e.  |  A/hot  I  =  IMcooil),  Pj\p  and  Pj{q  were  evaluated  to  be 
50  and  61.3%,  respectively,  of  P™ajX.  This  means  that  the  un¬ 
balanced  flow  rate  can  generate  larger  power  than  and 
P™c  shown  in  Fig.  5,  and  the  power  approaches  to  Pj™ 
when  the  flow  rate  of  a  fluid  increases  keeping  the  another 
constant.  The  fact  that  P^q  was  larger  than  P™p  is  consis¬ 
tent  with  the  flat  multi-panels  [11].  Therefore,  only  counter 
flow  type  was  considered  for  the  double  tube  systems. 

Pj2c h  was  96.0%  of  Pjiq,  while  T^ch  was  on^  35.2% 
of  Lpa£.  From  the  viewpoint  of  the  more  compact  system 
to  reduce  the  required  modules  and  the  space,  the  multiple 
tube  systems  are  more  effective,  and  the  T2CH  system  seems 
the  most  economic  for  construction.  These  tendencies  are 
similar  with  the  system  with  the  multiple  flat-panels  [11]. 
The  multi-tubes  systems  such  as  T3CH  or  T4CH  systems 
seem  worth  for  further  consideration. 

For  comparison  with  the  flat  panel  systems  in  the  pre¬ 
vious  study  [11],  Pmax  was  evaluated  in  Table  4  using  the 
identical  data  in  Tables  2  and  3.  Pmax  of  the  systems  with 
one  cylindrical  tube  (T1  types)  are  as  large  as  those  of  one 
flat  panel  systems  (FI  types).  However,  the  cylindrical  sys¬ 
tems  with  double-tubes  (T2CM  and  T2I)  cannot  generate 
the  larger  output  than  the  double  flat  panels  systems  (F2CM 
and  F2I),  except  for  the  helical  type.  The  reason  is  as  the  fol¬ 
lows.  The  over- all  heat  transfer  through  the  tube  is  roughly 
proportional  to  the  tube  radius.  Therefore,  the  heat  transfer 
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Table  4 


Output  power  and  system  length  when  the  output  power  becomes  maximum 


System  [11] 

Maximum  output 
power,  Pmax  (kW) 

System  length, 

Lmax  (m) 

System 

Maximum  output 
power,  Pmax  (kW) 

System  length 
Lmax  (m) 

F1P 

4.763 

18.72 

TIP 

4.744 

24.26 

F1C 

5.847 

29.80 

TIC 

5.825 

38.61 

F1I 

9.525 

37.44 

Til 

9.489 

48.51 

F2CM 

4.763 

37.44 

T2CM 

4.057 

35.74 

F2CH 

5.453 

12.39 

T2CH 

5.589 

13.64 

F2I 

8.594 

19.40 

T2I 

7.397 

14.29 

of  the  outer  tube  becomes  larger,  and  the  larger  amount  of 
heat  passes  through  the  outer  panel.  Therefore  the  fluid  tem¬ 
perature  changes  more  quickly  along  the  outer  path,  and  the 
temperature  difference  between  the  fluids  becomes  smaller 
at  the  outer  panel.  This  phenomenon  reduces  the  electromo¬ 
tive  force,  E,  especially  in  the  systems  such  as  T2CM  and 
T2I. 

^T2CH  was  2.50%  larger  than  as  listed  in  Table  4. 

This  is  because  the  helical  type  gives  a  homogeneous  tem¬ 
perature  difference  over  the  whole  panel  length.  This  homo¬ 
geneous  temperature  distribution  does  not  enhance  the  heat 
flow  through  the  outer  panel.  Then  the  temperature  differ¬ 
ence  on  the  thermoelectric  modules  become  larger  and  the 
output  power  increases.  The  helical  type  has  the  potential  of 
the  larger  output  by  increasing  the  number  of  panels. 

The  meandering  type  did  not  show  any  good  output  and 
miniaturization.  This  is  because  of  the  assumption  that  all 
the  thermoelectric  modules  were  connected  in  series  and 
that  all  panels  were  the  same  length.  If  we  can  control  the 
panel  length  of  two  panels  independently,  we  may  optimize 
the  output  power  of  the  individual  panel. 

6.  Conclusion 

This  work  studied  the  thermoelectric  power  generation 
with  cylindrical  multi-tubes  in  which  thermoelectric  ele¬ 
ments  were  embedded.  Six  systems  were  analyzed  using 
the  heat  transfer  theory  and  their  output  powers  were  ex¬ 
pressed  in  the  analytical  mathematic  equations.  Because  the 
over-all  heat  transfer  coefficient  depends  on  the  tube  ra¬ 
dius,  these  equations  are  more  complex  than  those  of  the  flat 
multi-panels. 

The  maximum  output  power,  Pmax,  was  obtained  at  the 
Til  system  where  the  cold  surface  was  kept  isothermally 
by  the  infinitively  large  heat  sink.  is  the  largest  in  the 
other  five  realistic  systems,  although  the  length  necessary  for 
the  maximum  output  is  longer  than  those  for  double  tubes 


systems.  The  helical  type  with  two  tubes  (T2CH)  can  gen¬ 
erate  as  large  as  the  TIC  system,  and  L^Ich  significantly 
shorter  than  that  for  TIC  system.  The  helical  type  with  the 
multi-tubes  makes  the  system  compacter,  keeping  high  Pmax. 
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